Measurements of sideward flow around the balance energy. 
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Sideward flow values have been determined with the INDRA multidetector for Ar+Ni, Ni+Ni 
and Xe+Sn systems studied at GANIL in the 30 to 100 A.MeV incident energy range. The balance 
energies found for Ar+Ni and Ni+Ni systems are in agreement with previous experimental results 
and theoretical calculations. Negative sideward flow values have been measured. The possible origins 
of such negative values are discussed. They could result from a more important contribution of 
evaporated particles with respect to the contribution of promptly emitted particles at mid-rapidity. 
But effects induced by the methods used to reconstruct the reaction plane cannot be totally excluded. 
Complete tests of these methods are presented and the origins of the "auto-correlation" effect have 
been traced back. For heavy fragments, the observed negative flow values seem to be mainly due to 
the reaction plane reconstruction methods. For light charged particles, these negative values could 
result from the dynamics of the collisions and from the reaction plane reconstruction methods as 
well. These effects have to be taken into account when comparisons with theoretical calculations 
are done. 



I. INTRODUCTION 



Studies of sideward flow, also called in-plane flow, have been found to provide information on the in-medium nucleon- 
nucleon interaction. By comparing the experimental results to dynamical calculations, it is possible to constrain the 
value of the in-medium nucleon-nucleon cross section a nn , and the incompressibility modulus of infinite nuclear 
matter Kqo The so-called balance energy Ebai (incident beam energy for which the sideward flow vanishes) 

has been found to be strongly dependent on a nn for light systems, and more dependent on Kqq for heavier systems 
0. For a fixed impact parameter and for a fixed incident energy, the flow parameter value strongly depends on 
Kqo HHQ. A dependence on the total isospin of the system has been also observed: keeping the total mass constant, 
higher values of Ebai are extracted for the more neutron rich systems ]To[ | . 

A simple interpretation of the sideward flow is that it results from first chance nucleon-nucleon collisions. At 
incident energies below Ebai , these collisions are sentitive to the attractive part of the nucleon-nucleon interaction 
and first chance particles are deflected toward an opposite direction relative to the initial projectile direction. In this 
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case, the flow parameter is negative. At energies higher than Ebai , the first chance collisions are sensitive to the 
repulsive part of the nucleon-nucleon interaction and the first chance particles are deflected toward the same direction 
relative to the initial projectile direction. In that case the flow parameter is positive. At Ebai , the repulsive and the 
attractive part of the interaction counterbalance and first chance particles are not deflected: the flow parameter is 
zero. 

But this scheme is indeed too simple. Dynamical calculations [fill have shown that the sideward flow may result 
from particles emitted at different stages of the reaction. Schematically, one contribution comes from the decay of the 
so called quasi-projectile and quasi-target, and the other one, from the emission of particles in the first moments of 
the collision. The relative rate between these two contributions is strongly dependent on the nature of the particle. 
This may explain the different values of flow actually observed for different types of particles P, p^ , ^3| . In this more 
realistic frame, the study of the detailed evolution of the sideward flow with the incident energy could shed a light on 
the production mechanism of particles around the nucleon-nucleon velocity. 

The measurement of small sideward flow parameter values (typically below 20 MeV/c/A) needs high accuracy and 
a complete information on each event which can be achieved by using powerful Air multi-detectors. Experimentally 
only positive values of the flow parameter can be measured, since the initial direction of the projectile is unknown, 
and since the positive direction is defined as the mean direction of particles emitted above the nucleon-nucleon frame 
velocity. A full understanding of the experimental methods is also needed, in order to correct possible spurious effects. 

The aim of this paper is to present the results of the sideward flow analyses on Ar + Ni, Ni + Ni and Xe + Sn 
systems from 25 to 95 A.MeV, and to make an extensive test of the standard methods used to measure the sideward 
flow. Sideward flow measurements for the Ni + Ni system at high energies can be found in E3. In the first section, 
the experimental set-up will be briefly described. The experimental results will be presented in the second section. 
The third section will be devoted to the test of various methods used to reconstruct the reaction plane. Conclusions 
will be drawn in the last section. 



II. EXPERIMENTAL SETUP 



The experiments were performed at the GANIL facility with the INDRA detector. Target thicknesses were respec- 
tively 193 /ig/cm 2 of 58 Ni for the 40 Ar + 58 Ni experiment, 179 ^g/cm 2 of 58 Ni for the 58 Ni + 58 Ni experiment and 
330 /ig/cm 2 of nat Sn for the 129 Xe + " at Sn experiment. Typical beam intensities were 3-4 xlO 7 pps. A minimal bias 
trigger was used: events were registered when at least three charged particle detectors fired. 

The INDRA detector can be schematically described as a set of 17 detection rings centered on the beam axis. In 
each ring the detection of charged products was provided with two or three detection layers. The most forward ring, 
2° < Oiab < 3°, is made of phoswich detectors (plastic scintillators NE102 + NE115). Between 3° and 45° eight rings 
are constituted by three detector layers: ionization chambers, silicon and ICs(Tl). Beyond 45°, the eight remaining 
rings are made of double layers: ionization chambers and ICs(Tl). For the Ar + Ni experiment the ionization chambers 
beyond 90° were not yet installed. The total number of detection cells is 336 and the overall geometrical efficiency 
of INDRA detector corresponds to 90% of Air. A complete technical description of the INDRA detector and of its 



electronics is given in 15 1(J ■ Isotopic separation was achieved up to Z=3-4 in the last layer (ICs(Tl)) over the 
whole angular range (3° < 9i a b < 176°). Charge identification was carried out up to Z=55 in the forward region 
(3° < Biab < 45°) and up to Z=20 in the backward region (9i a b > 45°). The energy resolution is about 5% for ICs(Tl) 
and ionization chambers and better than 2% for Silicon detectors. 

The INDRA detector capabilities allow one to carry out an event by event analysis and to determine reliable global 
variables related to the impact parameter. 

III. EXPERIMENTAL RESULTS 
A. Event sorting 

The first step was to sort events as a function of the violence of the collision. In this paper, we will use the total 
transverse energy: 

M 

E t rans,tot = H.lab X O"^)) 2 (1) 

i=l 
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where M is the charged particle multiplicity of the event, E\ lab , 9i are respectively the kinetic energy and the 
polar angle (with respect to the beam) of the particle i in the laboratory frame. QMD calculations have shown that 
Etrans,tot is a good indicator of the true impact parameter at intermediate energies fl7| . 

In order to sort events, we assumed a geometrical correspondence between E tra ns,tot and the impact parameter. 
The cross section of each bin was expressed as an experimentally estimated impact parameter, b exp jf8|. The first 
bin in E tra ns,tot, i.e the 31,4 mb of the highest values of E tra ns,tot, is linked to b exp between and 1 fm, the second, 
to b eX p between 1 and 2 fm, etc ... The last bin corresponds to b exp larger than 8 fm. This procedure has been 
applied on all detected events. Due to trigger conditions which remove the most peripheral reactions, the estimate 
of the impact parameter is not very accurate for the lowest values of E tran s,tot- From now, b eX r, will refer to the 
experimentally estimated impact parameter using the total transverse energy. Previous studies |J have shown that 
the flow parameter values are weakly sensitive on the exact choice of the sorting variable. 



B. Selection of "well measured" events 



The next step in the analysis was to select events in which sufficient information was recorded. This was achieved 
by requiring that the total measured J2i=i %i x Vpar (product of the charge Zj of particle i by its parallel velocity 
V par ) be larger than 70% of the initial Z pro j x V pro j of the projectile [Q. 

Since we want to study the dependence on impact parameter of in-plane flow, we checked that this selection 
conserves the whole impact parameter range. For all systems, we noted that the total transverse energy distribution 
of selected events covers the whole range of total transverse energy of registered events. If one assumes that the total 
transverse energy is a good measurement of the violence of the collision, then this result indicates that the whole 
impact parameter range of registered events is kept in selected events. Indeed, the major part of eliminated events 
corresponds to some peripheral collisions in which the target like fragment (TLF) was not detected, the projectile 
like fragment (PLF) was lost in the forward 0-2° beam hole and only few light particles were detected. In this case, 
Etrans.tot is still correctly measured and since the TLF and PLF transverse energies are very small. 



C. Flow parameters 



1. Definition 

To evaluate the flow parameters, one needs first to determine the reaction plane on a event by event basis. 
This is done by determining a transverse axis which defines with the beam axis the reaction plane. The transverse 
momentum method |l9| and the momentum tensor method pp| ] have been used. These methods have been found to 
be equivalent. As already mentioned in they give a better accuracy on the reaction plane determination than the 



azimuthal correlations method described in |21|. With such methods, the transverse axis is by definition oriented on 
the mean transverse direction of the forward emitted products (V par > V cm ). Within the standard interpretation of 
flow , the measured parameters values should be positive. More details about these methods will be given in section 
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In order to avoid the so called "auto-correlation" effect, the particle of interest is usually removed from the reaction 
plane determination, and a corrective momentum is added to the momentum of other particles [ p^| , p2|p^ ]. In this case, 
one plane per particle is determined. The tests of the methods and the issue of auto-correlations will be presented in 
section [v| 

Once the reaction plane is determined, the projection of transverse momenta on the reaction plane can be evaluated. 

For each reduced rapidity bin Y r = Y / Y pro j , its mean value <Pj; /A> is calculated, where Y — \log is 

the rapidity, v z the parallel to the beam component of the velocity in the laboratory frame and c the light velocity. 
The flow parameter F is by definition the increase of <P X /A> from Y=Y nn (Y r =0.5) to Y=Y pro j (Y r =1) by using 
the slope of the function <P X /A> =f(Y r ) at mid-rapidity |23| , p4]| . It reads: 

Z V ox T I (Y r =0.5) 



Typical evolutions of <P X /A> with Y r are shown in figures [j] and ||for Ar + Ni collisions at 74 A.MeV. The rows 
correspond to the particle types, the columns, to the experimental impact parameter bins. The momentum tensor 
method has been used to reconstruct the reaction plane. For each panel, open circles correspond to one plane per 
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particle (the particle of interest has been removed from the reaction plane reconstruction), full circles, to one plane 
per event (all products are taken into account for the reaction plane determination). For the Ar + Ni system, the 
balance energy is expected around 80 A.MeV for the central collisions. Therefore, small flow parameter values should 
be measured at 74 A.MeV. 

For the one plane per event prescription (full circles), the slopes at mid-rapidity are very high and much larger 
than the expected values. These overestimations are due to the so-called "auto-correlation effect" . The flow values 
determined with the one plane per particle prescription (open circles) are lower and, as a result, in a better agreement 
with the expected values. Some of them are even negative. Possible explanations of these negative slopes will be 
given in the following sections. 

2. Evolution of the flow parameters with the incident energy. 

For the most central collisions (h ex „ < 3 fm), the results are shown in figure || for the Ar + Ni system, in figure |] 
for the Ni + Ni system and in figure |3| for the Xe + Sn system. In the first two experiments, the reaction plane was 
determined by using the momentum tensor method with the one plane per particle prescription. The flow parameter 
values has been determined by a linear fit in the Y r range from 0.4 to 0.6. The errors on flow parameter values have 
been estimated by changing the Y r range from 0.35 to 0.65 and from 0.45 to 0.55. The error bar do not appear on 
the figures when it is smaller that the symbol size. 

The evolution of the flow parameter with the incident energy has a typical U shape for all particles. The incident 
energy which corresponds to the minimum flow value is located around 82 ± 2 A.MeV for the Ar + Ni system and 
around 75 ± 2 A.MeV for the Ni + Ni system. For both systems, these balance energies do not depend on the particle 
nature as for systems in reference [jTof| . They are in agreement with theoretical work 0. These calculation were 
performed with Koo w 220 MeV and a nn =0.8 c£^ ee , where cr£^ ee is the free nucleon-nucleon cross-section. 

For the Xe + Sn system (figure [5]), the results are shown for two reaction plane reconstruction procedures: the 
momentum tensor method (open diamonds) and the transverse momentum method (stars) with the one plane per 
particle prescription. As expected, both methods give close results. The U shape is barely seen for 3 He and heavy 
particles (Z > 3). For all particle types, the minimum flow energy is difficult to determine, since the flow parameter 
weakly depends on the incident energy. Unfortunately, the expected Eb a i value is around 50 A.MeV which is the 
maximum incident energy available for this system. No accurate determination of Eb a i can be done for this system. 
The full circles correspond to the momentum tensor method when one plane per event is determined. The flow 
parameter values are higher than the value obtained with one plane per particle due to the "auto-correlation" effect. 
The same behavior is observed with the lighter systems Ar + Ni and Ni + Ni when one plane per event is reconstructed. 

But the most striking feature of figures ||, || and || for the one plane per particle prescription is the observation of 
negative flow parameter values for d, t, 3 He and fragments with a charge greater than 3. By definition only positive 
values are expected. Negative flow parameter values have already been observed in previous studies [ p5| , but no 
clear explanation was given for this effect. We will propose in the next two sections two possible scenarios for these 
negative flow values. 

IV. THE PHYSICAL EFFECT FOR NEGATIVE FLOW VALUES. 

A possible explanation for the negative flow values is given by AMD calculations p~l]] : for light charged particles, 
the flow of promptly emitted (named "direct" ) particles can be opposite to the flow of "evaporated" particles (emitted 
from the quasi- projectile QP and the quasi-target QT). By definition, the reaction plane is oriented positively in 
the mean direction of forward emitted products. Therefore, negative flow values can be measured. 

The direction of the "direct" flow in the reaction plane results mainly from screening effects. As shown in figure 
, in the first moments of the reaction, "direct" nucleons coming from the projectile (target) are screened by the 
target (projectile) nucleus. For a fixed impact parameter, the orientation of the "direct" flow is determined by the 
geometrical configuration at the touching point and is weakly dependent on the incident energy. It is always aligned 
on the same direction as the initial orientation of the projectile in the reaction plane. At variance, the direction of 
the "evaporated" flow is aligned on the final directions of the QP. This direction depends strongly on the incident 
energy, as it will been explained in the next paragraph. Experimental studies have shown that heavy fragments 
could be emitted by non-evaporative processes, like a neck break-up |2^-|37]]. Since the direction of emission of these 
fragments is mainly aligned on the QP-QT axis, their flow direction is identical to the so-called "evaporated" flow. 
The contribution of heavy fragments to the flow parameter will be always attributed to the "evaporated" component. 
If one assumes that the observed flow parameter results only from these two contributions, the evaporated and prompt 
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emissions, a fairly simple picture can be proposed as shown in figures pL The simple pictures shown here are for 
pedagogical purpose. They do not pretend to reproduce the true process. They are shown to give a feeling of what 
kind of physical effect could lead to the observation of negative flow parameter values. 

At low incident energy, where the attractive part of the nucleon-nucleon interaction dominates, the QP and the 
QT are deflected to the opposite side relative to their initial directions (left column of figure ||). In this case, since 
<P X /A> is defined positive for particles which are deflected on the same side as the QP, the variation of <P X /A> with 
Y r is different for the two emissions: a "S" shape for the "evaporated" particles with a positive slope around mid- 
rapidity (dashed line in the lower left picture of figure ^) and a straight line with a negative slope for the "direct" 
particles (dashed and dotted line). 

At higher energies, where the repulsive part on the interaction is dominant, the QP "bounces" on the QT (right 
column of figure |^). In this case, the variations of <P X /A> with Y r give a positive slope around mid-rapidity for 
the two emissions as shown in the lower right panel of figure ^|. The observed variation of <P X /A> with Y r is a 
combination of these two contributions with their associated Y r distributions. 

Such combinations are schematically shown in figure 0. For each side panel, the original <P X /A> evolution of the 
"direct" particles (straight dashed and dotted line) and its corresponding Y r distribution (dashed and dotted gaussian 
distribution) are plotted. The <P X /A> contribution ("S" shape) and the Y r distribution of "evaporated" particles 
(sum of two gaussian distributions) are plotted as dashed lines. The resulting evolution of <Pa; /A> corresponds to 
the thick line. It is simply obtained by summing the two <Pa; /A> contributions weighted by their corresponding 
Y r distributions. 

If the contribution of "evaporated" particles is dominant, then the resulting flow parameter value is always positive 
(upper row of figure ^). The "S" shape is weakly affected by the contribution of "direct" emissions. At variance, 
if "direct" emissions become dominant (lower row of figure |7]) , then a negative flow value can be obtained. At low 
incident energies, the evolution of <P X /A> with Y r is dominated by the "evaporated" particles. When "direct" 
emissions becomes dominant, the evolution of <Pa; /A> with Y r follows the one of the "direct" particles, and the "S" 
shape is strongly deformed. Negative slopes can be found at mid-rapidity, i.e. negative flow parameter values. 

This explanation is very tempting for the light charged particles. Deuterons, tritons and 3 He are predominantly 
emitted in the mid-rapidity re gions, whereas protons and alpha particles are emitted as much by the QP and QT as 
by the mid-rapidity emissions |3qj39|| . At the same time, positive flow values are measured for alpha and protons and 
negative values for deuterons, tritons and 3 He (see figures |] and g). Both experimental features support the above 
interpretation. But this does not explain the negative values observed for the heavier fragments. Within the simple 
interpretation of sideward flow, the emission of fragments via cluster-cluster collisions of preexisting fragments in 
both partner bulk has a very low probability. As already mentioned, for PLF and the TLF, this observation is even 
in contradiction to the reaction plane orientation which is always oriented along the direction of the quasi-projectile. 
The solution of this puzzle has to be found elsewhere. 

Previous studies p3[ have shown that the reaction plane determination method and the associated issue of auto- 
correlations could strongly disturb the measurement of flow parameter. Before attributing the observation of negative 
flow parameter values to physical effects, one has first to be sure that these negative values are not due to the analysis 
methods. The next section is devoted to the test of reaction plane determination methods. The origins of the so 
called "auto-correlation" effects will also be studied. 



V. TEST OF THE REACTION PLANE DETERMINATION METHODS 

A. Reaction plane estimation. 

Let us now briefly describe three commonly used methods: the transverse momentum method [Q, the momentum 
tensor method |p(if and the azimuthal correlations method pl[ . From now on, the Oz axis is the beam axis, the Ox 
axis is the axis in the reaction plane which is perpendicular to the beam axis and Oy axis is the axis perpendicular 
to the reaction plane. The transverse momentum method, explained in details in |XS| ] , is based on the fact that the 
sum of transverse momenta of particles emitted by the quasi-projectile Pqp is opposite to the sum of the transverse 
momenta of particles emitted from the quasi-target Pqt- This is valid within the binary mechanism hypothesis, 
where the mid-rapidity contribution is negligible. Those vectors belong to the reaction p lane. I n ord er to maximize 
the efficiency of the method, one has to calculate the difference of the two vectors = Pqp — Pqt- Usually, (J is 
determined in the following way: 

N ^ 



5 



where N is the total number of particles in the event, Pj is the transverse momentum of particle j and uij a weight 
denned as follows: 

-1 if Y r - F rcm < -S 
~IY, ) = { if -5<Y X -y rcm <S (4) 
1 if Y x -Y Icm >5 



u{Y )=Y - Y rcm (5) 

where Y T is the reduced rapidity of the particle, Y Tcm the center of mass reduced rapidity (close to the reduced 
rapidity of the nucleon-nucleon frame Y rnn for the symmetric systems) and S a parameter which allows to remove 
the mid-rapidity particles from the estimation of the reaction plane. Q defines with the beam direction the reaction 
plane. With this definition, the reaction plane is systematically oriented along the quasi-projectile direction. 

Another way to estimate the reaction plane is to calculate and diagonalize a tensor T^ v . The beam axis defines 
with the eigenvector corresponding to the highest eigenvalue the reaction plane. As for the transverse momentum the 
reaction plane is oriented along the direction of the quasi-projectile. The tensor is defined in the following way: 

N 

= % X Pf x py with n,v = x,y,z (6) 
j'=i 

where P? is the momentum component along the fj, axis (/i = x, y, z) for the particle j. flj is a weight which is 
usually set to the inverse of the mass Aj of the particle flj — 1 /A,- (energy tensor) . 

The azimuthal correlations method is based on the following observation: in case of strong in-plane emission (high 
flow parameter value), the sum of the distances of particles momenta with respect to that plane are minimum 
This sum D 2 is calculated as follows: 



AT T / t->t , t->u\2 



^ 2 = E 



(7) 



where P* and P? are the transverse momentum components along the Ox and Oy axis respectively and a — tan(y>), 
where tp is the angle of the reaction plane relative to the Ox axis. Experimentally, one has to find the value of if which 
minimizes D 2 . With this method, the orientation of the reaction plane is not defined. One has to use the transverse 
momentum method to find it. In case of strong out-of-plane particle emission ( squeeze-out), the estimated reaction 
plane angle is wrong by 7r/2. 



B. Testing procedure 

Since the flow parameter is obtained from the average value of the transverse momentum projection on the reaction 
plane, one has to rebuild this plane from the experimental data. We have checked the reliabilities of the three reaction 
plane reconstruction methods. The general procedure of the test is the following: a known flow parameter value F is 
set for a sample of generated events; then the reaction plane estimation method is applied on that sample and the so 
called "experimental flow parameter" F exp is determined. A method is considered effective if the experimental value 
F exp is equal or close to the initial one Fi n u ■ This allows to also check the additional disturbance introduced by 
the experimental set-up compared to the method itself. 

To set the flow parameter value F to an event, an in-plane component P£ is added to the transverse momentum 
of each particle, similarly to the procedure used in p3|. The amplitude of this in-plane component depends on the 
reduced rapidity of the particle: 

if y r < 0.25 
0.5) if 0.25 < Y T < 0.75 (8) 
if Y t > 0.75 
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where F is the flow parameter value to be set on, Y r is the reduced rapidity of the particle, A its mass and P£(F, Y r ) 
the in-plane component. 

The test has been performed on two systems Ar + Ni at 74 A.MeV and Xe + Sn at 50 A.MeV. These two systems 
have been studied with INDRA and their comparison allows to check the effect of the mass of the system on the 
transverse flow measurement. These energies values have been chosen because they are close to the expected balance 
energy. For each system, 15000 events have been generated using the SIMON code, whose entrance channel includes 
a pre-equilibrium emission of protons and neutrons JfCj ] . SIMON is not used to reproduce the experimental data, 
but rather as an event generator. The goal is to check how the reaction plane reconstruction methods react in a well 
defined situation. In this test, only the most central collisions have been used. The flow parameters for the different 
particle types are not exactly zero but close to zero, and differ from a particle type to another. The effect of the 
addition of the P£ component is to add the value of F to the original value of the flow parameter, giving a F init 
value. This procedure will allow to study the effects of the different reaction plane reconstruction methods on flow 
parameters measurements, as in [ p3| . 



C. One plane per event or one plane per particle? 



Since the transverse momentum is used both for the reaction plane estimation and for the projection, auto- 
correlations are expected (see full circles in figures [j] and ||). Auto-correlation effects are also amplified by the 
loss of information due to a non perfect detection. The usual way to solve this problem is to remove the particle 
which has to be projected from the estimation of the reaction plane. Thus, an additional corrective component pj°^ 
is added to the momentum of the remaining particles in order to ensure the conservation of the total momentum. Its 
definition is the following: 



pco 

j 



f 



A; 



2^fc=l 



A L 



V, 



(9) 



where i is the removed particle, Vi its velocity, Aj the mass of particle j (j ^ i) and M the event multiplicity. In 
this case, one plane is determined for each particle of the event. 

Two prescriptions have been used for the methods described in section |V A| . In the first one, one plane per event 
is determined. This allows to check the effect of the auto-correlations with respect to the method used. In the second 
one, one plane per particle is determined, in order to test the efficiency of the correction. 



1, Azimuthal dispersion between the true reaction plane and the reconstructed one 

To compare the relative efficiencies of these methods, the distribution of the angular azimuthal difference A$ 
between the true and the reconstructed reaction plane directions has been studied. Such distributions have been 
already shown (see figure 6 of §|). 

The observed mean value < A$ > is zero. The accuracy of the reaction plane determination is estimated with 
the standard deviation u(A$) of these distributions as a function of the added flow parameter value F. Figure || 
shows such evolutions of cr(A$) with F. F has been used instead of Fi n a because the Fi n n values are different from a 
particle type to another, whereas the same F value has been added to all particles. The upper row corresponds to the 
Ar + Ni system, the lower row, to the Xe + Sn system. The left column correspond to the case when one plane per 
particle is determined and the right column to the one plane per event prescription. No significant difference is found 
between these two prescriptions. For all methods and systems, the dispersion tT(A$) decreases when F increases. For 
small F values, the different methods give slightly different results because the initial flow parameter is not zero. As 
expected, the reaction plane determination is more accurate in case of strong in-plane emission. For both systems, the 
transverse momentum method and the tensor method give similar accuracies, whereas for the azimuthal correlation 
method the dispersion is systematically higher for all F values. For the Xe + Sn system, the values of <ta$ are smaller 
than for the Ar + Ni system. 

From these observations, three main conclusions can be drawn: i) the azimuthal dispersion of the reconstructed 
reaction plane is the same for the one plane per particle and one plane per event prescriptions; ii) similarly to the 
conclusion of reference ](|, the azimuthal correlation method is less accurate than the two other methods, even at low 
Finn values; iii) the higher the mass of the system and the higher the Fi n a value, the more accurate is the reaction 
plane reconstruction. 
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Surprisingly, no difference is seen on cr(A$) between the one plane per particle and the one plane per event 
prescriptions, whereas strong discrepancies are seen for the <P X /A> =f(Y r ) curves (figures |l| and One has 
to keep in mind that in the present case, only the deviation from the true reaction plane is studied. For the 
<P X /A> =f(Y r ) curves, the combined effects of the reaction plane reconstruction accuracy and of the projection of 
transverse momenta on the reconstructed reaction plane are present. 

2. Measured flow parameter versus initial flow parameter 

For the transverse momentum method, the dependence of F exp on the initial value of the flow parameter is shown 
in figure ^| in the case of one plane per particle (left column) and in case of one plane per event (right column) . For 
one plane per particle, the procedure used to remove auto-correlations seems to be efficient for H isotopes. For 
heavier particles, F exp is systematically below Fi n n . This underestimation increases with increasing charge. For a 
given Fi n it value, the underestimation is lower for Xe+Sn system than for Ar+Ni. Finally, the amplitude of this 
underestimation decreases with increasing value of F^a . 

With one plane per event (right column of figure ||), F exp values are systematically larger than Fi n u values due 
to the auto-correlation effects mentioned above. For a fixed Fi n u value, the overestimation increases with decreasing 
charge of the particle. Here again, the amplitude of the overestimation diminishes for higher values of Fi n u and with 
increasing mass of the system. 

For the momentum tensor method, results are shown in figure ^ for one plane per particle (left column) and for 
one plane per event (right column). The same trends are observed as for the transverse momentum method. The 
auto-correlation effects are smaller in the case of one plane per event. 

For the azimuthal correlation method (figure [ll]), the same trends are observed as for the two other methods. The 
main difference is a somewhat larger underestimation of the flow parameter value. 

3. Influence of the experimental set-up 

Part of the anti or auto correlation may be related to the detector efficiency. In order to probe the effect of the 
experimental set-up, the INDRA "filter" has been applied on the generated events. The INDRA "filter" is a 
sophisticated software which simulates the response of the detector. For each particle, the energy losses in each layer 
of the detector are calculated. An identification procedure similar to the experimental one is then applied on the 
energy losses giving back the charge and the energy of the detected particle. Doing this, multiple hits in a detector 
are treated in the same way they are treated in the experiment. This procedure allowed to reproduce the angular and 
energy thresholds observed experimentally. 

The reaction plane is then reconstructed from the so called "filtered" events. In figure |l2|, the correlation between 
F exp and F% n it is plotted for the momentum tensor method with the one plane per event prescription. For the two 
systems studied here, the results are very similar to those obtained with a perfect detection (see the right column of 
figure [n^) . The measured flow values are above the initial ones for the Ar + Ni system, and close to the initial ones 
for the Xe + Sn system. The detector has a weaker effect than the reaction plane determination procedure. This 
conclusion is identical to those made in reference p^] for a 4ir array which had higher thresholds. 

4- Conclusions of the simulations 

The following conclusions can be drawn from the previous study: 

i) removing the particle of interest from the reaction plane introduces an anti-correlation which is not counterbalanced 
by adding a corrective momentum Pj°^ to the momentum of other particles. The anti-correlation effect leads to an 
underestimation of sideward flow values; 

ii) when one plane per particle is determined, the amplitude of the anti-correlation is higher for small Fi n u values, 
for heavy particles and for small system size; 

iii) when one plane per event is determined, the amplitude of the auto-correlation is higher for small Fi n n values, for 
light particles and for small system size; 

iv) for all methods, for the one plane per event prescription, the auto-correlation effect leads to an overestimation of 
the flow parameter value. This overestimation is minimum for the momentum tensor method; 

v) the best results are obtained for the momentum tensor method with the one plane per event prescription; 

vi) Detection effects are weaker than effects induced by the procedures used to reconstruct the reaction plane. 
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One could conclude from these simulations that the best method to measure small sideward flow values is the 
momentum tensor method with one plane per event. Unfortunately, the simulation is too simple compared to the 
experimental situation. We remind the reader that in this simulation, the mid-rapidity contribution is only present 
for protons and neutrons. This is why the "auto-correlation" effect is mainly seen for protons in the simulation. In 
the experimental data, the mid-rapidity emission is also made of heavier particles |3§|]39|] . The autocorrelation effect 
is seen for all particles for the Ar + Ni system (full circles in figures [l] and ||) as well as for the Xe + Sn system (full 
circles in figure ^|). In the latter case, the flow parameter is weakly dependent on the incident energy. Therefore, 
from the experimental results, it turns out that the "auto-correlation" effects are strong. In other words, qualitative 
effects may be understood thanks to the simulations but quantitative estimations of the corrections required in the 
data are difficult to realize from these simulations. Such a correction may be done if the "auto-correlation" effect is 
well understood. The study of the origins of the "auto-correlation" effect is done in the next section. 



D. Origin of the auto-correlations 



These studies show us that the methods developed at high energies, where high values of sideward flow are 
measured, are not well suited for intermediate energies where the flow parameter values are typically around or below 
30 MeV/c/A. The amplitude of the anti or auto correlations depends also on the nature of the particle. Let us try 
nevertheless to identify the origin of the auto-correlations at intermediate energies. 

First of all, one has to make some remarks. If a method would be able to reconstruct perfectly the reaction plane 
from the momentum of all particles, no auto-correlation effect would be seen. This is shown for example in the right 
column of figure [l^ for the Xe+Sn system. The experimental flow parameter F exp is very close to Fi n u for Fi n u values 
above 20 MeV/c/A, although the momenta of all particles have been used to calculate the momentum tensor. This 
is quite unexpected since for high Fi n u values, the transverse momenta values are large. The auto-correlation effect 
should be maximum for high Fi n u . The usual explanation of the auto-correlation effect does not seem to be the 
right one. If so, where does this "auto-correlation" effect come from ? Since the azimuthal correlation method is less 
accurate than the two other ones, the origin of the "auto-corellations" will not be checked for this method. 



1. "Auto- correlations" for the transverse momentum method 



In the transverse momentum method, one assumes that the particles emitted above Y cm are all coming from the 
decay of the quasi-projectile (QP), and those emitted below Y cm are all coming from the decay of the quasi-target 
(QT). But at intermediate energies, the contributions from the QT, the QP and from the mid-rapidity area are 
mixed, especially for the most violent collisions |3^,|39]]. In addition, for the most central collisions, the azimuthal 
angular distributions are rather flat and no privileged direction can be clearly seen. Therefore a wrong weight may 
be attributed to the particles, and the estimated reaction plane may have nothing to do with the true one. In this 
case, the reconstructed reaction plane may be oriented along the particles with the highest momenta. 

On the other hand, if the right weight was attributed to the right particles, this effect should vanish. This can 



be checked in the simulation, for which the origin of particles is known. The results are shown in figure 13, In 



this simulation, a pure binary scenario has been assumed: the first stage of the collision leads to the formation of a 
quasi-projectile and a quasi-target both deflected in the reaction plane, without any pre-equilibrium emissions. One 
plane per event is determined, using the transverse momentum method, and the weight of the particles is attributed 
according to their origin: +1 for the particles emitted by the quasi-projectile, -1 for those emitted by the quasi-target. 
It is seen that the effect of "auto-correlation" is removed, even for the smaller flow values. The so called "auto- 
correlation" effect comes from a loss of information (the origin of the detected particles), instead of the use of the 
transverse momenta in both the reaction plane determination and in the projection. 

In the experiment, the exact knowledge of the origin of the particle is impossible, especially for the most damped 
reactions. In addition the collisions are not purely binary due to prompt emissions. The promptly emitted particles 
carry a part of the total tansverse momentum and the QP and QT are hence pushed out of the reaction plane. In the 
experiment, the perfect determination of the reaction plane using the tranverse momentum method is very difficult, 
especially for the central collisions. 
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2. "Auto- correlations" for the momentum tensor method 

For the momentum tensor method, the origin of "auto-correlations" effects can be understood using a simple test. 
Let us consider the case where a quasi-projectile of mass number A splits into two equal size fragments and the 
quasi-target of the same mass number A remains unchanged. The quasi-projectile and the quasi-target are deflected 
in the reaction plane. The axis joining their center of mass has an angle Qdefiec with respect to the beam axis (see 
figure The axis joining the two fragments issued from the splitting of the quasi-projectile has an angle 9 sp ut with 
respect to the beam direction and an angle 4> sp ut with respect to the reaction plane. The relative velocity between the 
quasi-projectile and the quasi-target is V r and the relative velocity between the two fragments of the quasi-projectile 
is aV r . A scheme of the described configuration is presented in figure [l4| 

The momentum tensor of such a simple case can be calculated and one can study the azimuthal angular difference 
A0 between the reconstructed reaction plane and the true one. The direction of the reaction plane is defined in the 
plane transverse to the beam axis, displayed in the right panel of figure [l4| In this panel, the QT is by definition in 
the reaction plane but the two splitting fragments are out of the plane. 

The first key variable is a. When its value is small (close to zero), the QP-QT axis is obviously the main axis and 
the reaction plane is perfectly determined. But when the a value is big enough the axis between the two splitting 
fragments may become the main axis. In this case, the reconstructed reaction plane direction may be dependent on 
the splitting direction (f> S pUt which is the second key variable. 

In figure |l5|, the evolutions of Ac/) as a function of a and <fi S piit are displayed. The angles Qdefiec and 9 sp ut are 
set respectively to 10° and 40° . Similar pictures are obtained for other 9d e fiec and 9 sp ut values. For large enough 
values of a A0 is strongly dependent on 4> S pUt- For low values of a, i.e. small relative velocities of the two fragments 
compared to the relative velocities of the quasi-projectile and the quasi-target, A<f> does not depend on <fi S piit and is 
equal to zero. The two out-of-plane fragments do not introduce much pertubations on the reconstruction procedure. 
The direction of the reconstructed reaction plane is mainly determined by the QT. The reaction plane is therefore well 
estimated. On the other hand, when a is larger than 1, the reaction plane is mainly determined by the two out-plane 
fragments. The direction of the reconstructed reaction plane is therefore correlated to the splitting direction. In the 
experiment, this last configuration is similar to the most violent (central) collisions but with a larger multiplicity. 

If the two fragments issued from the splitting of the quasi-projectile were gathered before applying the momentum 
tensor method, the reaction plane would be perfectly determined whatever the fragmenting configuration. That 
means that if the origin of the fragments is known, on can determine perfectly the reaction plane by grouping the 
fragments coming from the same source in a single fragment. The results of the momentum tensor method depend 
on the way the fragments are gathered. As for the tranverse momentum method, a perfect determination of the 
reaction plane by using the momentum tensor method is very difficult in the experiment. 

E. Discussion 

For incident energies around the balance energy, the standard methods used for the reaction plane reconstruction 
are not well suited. Whatever the method used, the one plane per event prescription leads to an "auto-correlation" 
effect i.e. a large overestimation of the flow parameter values. At variance, the one plane per particle prescription 
induces an anti-correlation effect which gives an underestimation of the flow parameter. 

This indicates that the negative flow values observed in experimental data can be attributed to the used recon- 
struction methods, especially for the heaviest fragments. For the light charged particles, the physical effect can not 
be completely ruled out, since positive values of flow are observed for alphas, whereas negative values are expected if 
the reaction plane determination methods effects are dominant. For such particles, the two effects are probably mixed 
and more detailed studies have to be performed to establish their relative weights in the observed values. To obtain 
the true flow parameter values, one has to understand the auto-correlation effects in order to correct them accurately. 

But understanding the origin of the "auto-correlation" effects is a complicated task. To correct them, a complete 
knowledge of the origin of particles is needed. This can be achieved only for the less violent collisions and/or at 
higher energies where the mixing between the different contributions is weak. For the most violent collisions, such a 
knowledge is unreachable unless assumptions are made. But in this case, the flow values obtained may only result 
from these assumptions. 

Since the correction of experimental data seems to be impossible, it could be easier to apply the experimental filter 
and the analysis procedure on theoretical calculations. Most of the available dynamical calculations have to evolve to 
enable this procedure, since most of them are following the time evolution of the one body density. More precisely, 
the dynamical calculations should include the proper description of particle and fragment formation. 
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VI. CONCLUSIONS 



The in-plane flow parameter has been determined for the Ar + Ni collisions from 32 A.Mev to 95 A.Mev, for the 
Ni + Ni system from 32 A.Mev to 90 A.MeV and for the Xc + Sn system from 25 A.MeV to 50 A.MeV. For central 
collisions, the balance energies are equal to 82 ± 2 A.MeV for the Ar + Ni system and 75 ± 2 A.MeV for the Ni + Ni 
system. For the Xe + Sn system, the balance energy is around 50 A.MeV, but experiments at higher incident energies 
have to be performed to confirm the result. These values are in agreement with the systematics of balance energies 
found in other experiments. This systematics has been reproduced by a dynamical model assuming Koo w 220 MeV 
and a nn =0.8 cr^ ee @]. As already observed, the balance energy weakly depends on the particle nature. For these 
central collisions, negative flow values are observed for both systems. Two different explanations are proposed for 
these negative values. 

The first one, supported by transport model calculations, attributes this effect to the relative importance between 
the prompt emission and the evaporative one. A negative flow parameter value can be observed if the prompt 
emission is dominant. This explanation seems to be satisfactory for the light particles. Deuterons, tritons and 3 He 
are predominantly emitted in the mid-rapidity regions []38| , j39| and their flow parameter values are indeed negative. 
At variance, for protons and alpha particles, the measured flow parameters are positive. They are emitted as much 
by the QP and QT as at mid-rapidity. But on the other hand, negative flow values are measured for fragments 
while the prompt emission is not the dominant process for these products. For these heavy fragments, the observed 
negative values cannot be explained in this way. 

The second explanation attributes these negative values to the experimental methods used to extract the reaction 
plane. The usual method used to avoid auto-correlations, the omission of the particle of interest, leads to an anti- 
correlation. This induces an inversion of the reaction plane reaction and then lead to the measurement of negative 
flow values. The amplitude of this effect increases when the flow parameter value decreases i.e. when the incident 
energy is getting closer to the balance energy. This explanation is supported by the observation of negative flow values 
for the heaviest fragments, whereas a positive value is expected. A careful study of the "auto-correlation" effect shows 
that its manifestation results from the loss of information about the product origins for both methods. 

In experimental data, these two effects are probably mixed up. They disturb the measurement of the absolute 
value of sideward flow, especially around the balance energy for which low flow parameter values are expected. On the 
other hand, the relative evolutions with incident energy, and especially the determination of the balance energy, are 
in agreement with previous experimental studies and theoretical calculations. It may indicate a relative robustness 
of the balance energy variable. In the present status, the real effect can only be studied with simulations on which 
the complete experimental procedure can be applied. An accurate determination of the in-medium nucleon-nucleon 
interaction parameters can only be achieved if the disturbances induced by the analysis methods and the experimental 
set-up are explicitely taken into account. This requires an evolution of dynamical calculation to make possible this 
comparison procedure. 
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FIG. 11. Same as Blby using the azimuthal correlation method. 
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FIG. 12. Same as RO right column but the events have been affected by the INDRA filter. 
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FIG. 13. Correlation ol the measured flow parameter value F exp and the initial flow parameter value Finit in the framework 
of the SIMON code assuming a pure binary scenario (excited QP and QT only). Left column: the transverse momentum 
method has been used and the weights are determined according to the reduced rapidity of the particle. Right column: the 
weights are set according to the true origin of the particle. One plane per event has been calculated 
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